Ising Deconfinement Transition Between Feshbach-Resonant Superfluids 
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We investigate the phase diagram of bosons interacting via Feshbach-resonant pairing interactions 
in a one-dimensional lattice. Using large scale density matrix renormalization group and field 
theory techniques we explore the atomic and molecular correlations in this low-dimensional setting. 
We provide compelling evidence for an Ising deconfinement transition occurring between distinct 
superfluids and extract the Ising order parameter and correlation length of this unusual superfiuid 
transition. This is supported by results for the entanglement entropy which reveal both the location 
of the transition and critical Ising degrees of freedom on the phase boundary. 
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The ability to cool atoms to low temperatures, and 
control their interactions, has revolutionized the study 
of quantum many body systems. Important achieve- 
ments include realizations of Bose-Einstein condensa- 
tion (BEC), Bardeen-Cooper-Schrieffer (BCS) pairing 
in Fermi gases, and strongly correlated Mott insulators 
(Mis). In this development, the BEC-BCS crossover be- 
tween a gas of tightly bound molecules and weakly bound 
Cooper pairs has played an instrumental role, and it has 
been widely explored using Fcshbach resonances to in- 
duce pairing. This has led to diverse studies of the 
condensate fraction, single particle gap, collective exci- 
tations, and vortices, and to pioneering approaches to 
molecular quantum chemistry. For a review see Ref . [l| . 

In recent work it has been argued that the BEC- 
BCS "crossover" for bosons is strikingly different to 
the fermionic case since the atoms as well as molecules 
may undergo Bose-Einstein condensation. These stud- 
ies have raised the exciting possibility of an Ising quan- 
tum phase transition between distinct molecular (MC) 
and atomic plus molecular (AC+MC) condensates. In 
addition to discrete Z 2 symmetry breaking, this transi- 
tion has a topological character and may be viewed as a 
confinemcnt-dcconfincmcnt transition for vortices. 

The principal aim of this manuscript is to establish the 
presence of such novel Z2 transitions in one-dimensional 
(ID) bosonic Feshbach systems, where strong quantum 
fluctuations destabilize long range superfiuid order. We 
combine large scale density matrix renormalization group 
(DMRG) Q and field theory techniques to provide com- 
pelling evidence for Ising behavior. We elucidate a full 
characterization of the scaling regime and the proximate 
phases. Our results demonstrate that an Ising transi- 
tion survives at strong coupling and large densities where 
field theory arguments are no longer justified. For re- 
lated transitions in the attractive Bose-Hubbard model 
with three-body losses see Refs. [1, El], and for analogues 



involving multicomponent fermions see 1^, U 
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FIG. 1. Phase diagram of the ID Hamiltonian |T]) with total 
density pt = Nt/L = 2, showing a Mott insulator (MI), 
a molecular condensate (MC), and a coupled atomic plus 
molecular condensate (AC+MC). We use DMRG with up 
to L — 128 and open boundaries, with e a = 0, Uaa/2 = 
U m m/2 = U am = g = U, t a = 1, t m = 1/2. The squares 
and circles indicate the vanishing of the one-particle and 
two-particle gaps, E\ g and E 2g , as L — > oo. The crosses 
show where the molecular exponent, i/ m , reaches unity. Inset: 
AC+MC to MI transition at e m — 4. The atomic and molec- 
ular exponents, v a and v m , are locked up to the MI boundary 
where v m — 1, indicating the absence of an AC phase. 



We consider the Hamiltonian 0, @, [H| 



r . Tli a Tli 



(1) 



describing bosons, bi a , hopping on a lattice with sites i, 
where a = a,m labels atoms and molecules. Here, e a are 



2 



on-site potentials, t a are hopping parameters, (ij) de- 
notes summation over nearest neighbor bonds, and U aa ' 
are interactions. Normal ordering yields : n ia n ia := 
nia( n ia ~ 1) f° r like species, and : n ia n ia r := n ia n ia i 
for distinct species. Molecules are formed by the Fesh- 
bach term, Hp = gJ2i( m l a i a i + H.c), where m, = 6j m 
and a,i = bia- Atoms and molecules are not conserved, 
but the total, iVr = J2i( n ia + 2ni m ), is preserved. 

To orient the discussion, we present a section of the 
phase diagram in Fig. [TJ with parameters chosen for com- 
parison with previous studies @. In this manuscript we 
use DMRG on a ID system with up to L = 512 sites, and 
adopt units where t a = 1. We allow up to five atoms and 
five molecules per site, and retain up to m p = 2400 states 
in the density matrix so that the discarded weight is less 
than 1 x 10~ 10 . The phase boundaries correspond to the 
vanishing of the one-particle and two-particle excitation 
gaps, Ei g ee (ii p (L) - fiih(L) and E 2g = P2 P {L) - (J*2h(L) 
respectively, where p np {L) = E (L, Ar + n) — i? (£, iVx), 
/j, nh (L) = E Q (L,N T ) - E Q {L,N T - n), and E Q is the 
ground state energy. The diagram shows a MI with gaps 
for both excitations E\ g + 1 and E 2g ^ 0, a MC phase 
with a one-particle gap E\ g ^ and E^g = 0, and a cou- 
pled atomic plus molecular condensate (AC+MC) with 
E\ g = and Eon =0. In contrast to the qualitative 
diagram in Rcf. [6( , inferred from quantum Monte Carlo 
simulations on smaller systems, we find no evidence for 
a single-component AC phase. This is in accord with 
expectations in higher dimensions [2[. As we will dis- 
cuss, this is supported by direct evaluation of correlation 
functions using both DMRG and field theory. Through- 
out the AC+MC phase we find power laws for atoms and 
molecules with related exponents; see inset of Fig.[TJ The 
conclusions of Ref. Q are hampered by the slow diver- 
gence of the associated zero momentum molecular occu- 
pation number with increasing L, close to the MI bound- 
ary. This also afflicts the molecular visibility. Here, our 
focus is on the transition between the MC and AC+MC 
superfluids. We begin with symmetry arguments and 
field theory predictions before comparison with DMRG. 

An intuitive way to understand the origin of the pro- 
posed Ising transition between the MC and AC+MC 
phases is via the symmetry of the Hamiltonian ([T]) un- 
der U(l) x Z 2 transformations. This corresponds to in- 



variance under to 



to and a 



f (0/2±7T) a wnere 



9 € R. In general these symmetries may be broken inde- 
pendently. Before discussing the problem in ID, where 
continuous U(l) symmetry breaking is absent, let us first 
recall the situation in higher dimensions Q • In this case, 
the molecular condensate (MC) phase has (m) ^ and 
(a) = 0. This only breaks the U(l) contribution and 
leaves the Z 2 symmetry, a — > —a, intact; this corresponds 
to the disordered phase of an Ising model, coexisting with 
molecular superfluidity. On the other hand, the coupled 
atomic plus molecular condensate (AC+MC) phase has 
(to) 7^ and (a) + 1 0. This breaks the U(l) x Z 2 symme- 



try completely and corresponds to the ordered phase of 
an Ising model, coexisting with atomic and molecular su- 
perfluidity. Returning to the present ID problem, where 
continuous U(l) symmetry breaking is absent, the spon- 
taneous formation of expectation values (a) and (m) is 
prohibited. Instead, supcrfluid order is characterized by 
long range power law correlations, and the nature of the 
phases and transitions in Fig.[T]requires closer inspection. 

Owing to the U(l) x Z 2 symmetry of the Hamiltonian, 
the low energy Lagrangian of the MC to AC+MC tran- 
sition is given by C = Cg + £^ + C§^> 0, HJ , where 



[c^ 2 {d T df + {d x vf 



is a U(l) invariant free scalar field, and 



r _K 4> 



c^\d T <t>) 2 + {d x <t>) 2 -r,cj> 2 + \cf> 



(2) 



(3) 



is an Ising model in the soft-spin </> 4 representation. The 
coupling, Lfltj, = icj) 2 d T '&/2, has a similar form to a Berry 
phase [2; 3- A similar action also emerges for quan- 
tum wires [l3|. In the following we neglect L^^ and 
examine the reduced theory. Within mean field theory, 
£i9</> ~ i((j)) 2 d T d /2 acts like a boundary term, and this is 
expected to provide a good description of the proximate 
phases. Near the transition, this cannot be neglected a 
priori, and Ls^ may change the behavior on very large 
leng th scales and in other regions of the phase diagram 
|l3[. Nonetheless, we find excellent agreement with bulk 
properties. The parameters K$, c$, K<p, c,p, rj, A, are re- 
lated to the coefficients of H . Atoms and molecules are 
described by the semiclassical number-phase relations, 
m ~ \J~Pm e rf , and a ~ <j> e lf> ^ 2 , where p m is the molecular 
density. We will explore the consequences of this corre- 
spondence in ID, for local observables and correlations. 

Let us first gather consequences of this correspon- 
dence for local observables. Deep within the Z 2 dis- 
ordered MC phase, rj ^> and ((f>(x)) = 0. How- 
ever, 4> 2 { x ) rnay have a non-zero average. It follows 
that the densities of atoms and molecules, (a^ (x)a(x)) ~ 
((j) 2 (x)) and (m^ (x)m(x)) ~ p m , are generically non- 
zero in both the AC+MC and MC phases. In addi- 
tion, (m) (x)a(x)a(x)) ~ yjp m ( <fi 2 (x)), acquires true long 
range order, even in this ID setting; Hp locks the atomic 
and molecular condensates as encoded in the number- 
phase relations. However, this local average is naively 
insensitive to the Z 2 transition due to invariance under 
a —> —a. Insight is better gleaned from correlations. 

that 



ii) 



It follows from the relation to ~ \[Pm e 
the molecular correlation function (rrv (x)m(0)) ~ 
Pra(e _! ''' I 'e"'' ') ~ x~ Vm decays like a power law, 
where v m = 1/2ttK^ varies throughout the phase dia- 
gram. In contrast, the behavior of the atomic correlation 
function, (ot(a;)o(0)} ~ (<f>(x)<f>(0)) ( e -«(*)/2 e «(°)/ 2 ) ~ 
((/>(x)(/)(0))a; _l/m / 4 , depends on the Ising prefactor. We 
consider the disordered and ordered phases in turn. 
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FIG. 2. Correlation functions using DMRG with L = 128 
and open boundaries. We use the parameters in Fig. [T]with 
U = 0.7. (a) Z2 disordered MC phase with e m — —4, re- 
vealing power laws for molecules and atomic bilinears with 
the same exponent; the fits are y = 0.858 x~ 0,1922 and 
y = 0.130 a; -0,1909 . Inset: Atomic correlations decay exponen- 
tially. We fit to the prediction (a t (x)a(0)} ~ x~ Um/4 K {x/£), 
where we input v m from panel (a) and extract £ ~ 9.28. This 
establishes MC as a pairing phase without atomic condensa- 
tion, (b) Z2 ordered AC+MC phase with e m = — 3. Atoms 
and molecules exhibit power law exponents locked by a factor 
of four; the fits are y = 0.667 a;" 1827 and y = 0.657 z~ ' 0456 . 



In the Z2 disordered MC phase, the atomic correlation 
function decays exponentially with a power law prefac- 
tor, (at(x)a(x)) ~ i-"»/*Ko(x/0 ~ s -i/a-«'m/* e -*/C. 
Here we use the result for the hard-spin Ising model, 
((j>(x)<fr(0)) ~ K (x/£), where K is a modified Bessel 
function and £ is the Ising correlation length [l4[ . On the 
other hand, pairs of atoms condense and exhibit power 
law correlations, (af (x)a^ (x)a(0)a(0)) ~ (0 2 ) 2 x~ Vm , 
with the same exponent as the molecular two-point func- 
tion, v m . That is to say, the MC phase is a pairing phase 
of bosons without single particle condensation [151 ] . In 
order to test these weak coupling predictions we per- 
form DMRG on the ID Hamiltonian (JXJ) . As predicted, 
this behavior is well supported by our simulations in 
Fig.[5fa), which reveal identical power laws for molecules 
and atomic bilinears, with exponential decay for atoms. 
This behavior extends throughout the MC phase, includ- 
ing the Mott boundary in the strongly coupled regime. 

In contrast, in the Z2 ordered AC+MC phase, 
both molecules and atoms have power law correlations, 
(mt(x)m(O)) ~ a;""-, ( O f(s)o(0)) ~ (<f>) 2 x-^, with 
locked exponents, v m = 4v ai that differ by a factor of 
four Again, these features are readily seen from our 

large scale DMRG simulations in Fig. H^b) and the inset 
of Fig. [T] Likewise, this behavior persists into the strong 
coupling limit, where the field theory approach no longer 



strictly applies. In particular, we have checked that the 
molecular correlation function, (m'(x)m(0)) ~ x~ Vm , re- 
mains a power law throughout the AC+MC phase and 
close to the Mott boundary in Fig. [T] This is consis- 
tent with the absence of an AC phase Q in contrast to 
Ref. 0. The latter employ the zero- momentum occupa- 
tion, n(0). However, the Fourier transform of x~ v gives 
n(0) ~ const + const L 1 ~ v \ close to the MI where v m = 1 
one may miss the slow divergence of n m (0). 

Having established a close connection between field 
theory and DMRG for the MC and AC+MC phases, 
let us now examine the transition. A key diagnostic 
is the central charge, c, which counts critical degrees 
of freedom. This may obtained from the entanglement 
entropy For a block of length I in a periodic sys- 
tem of length L, the von Neumann entropy is given by 
Sjr(Z) = —Tti(pilnpi), where pi = Tr L _ ; (p) is the re- 
duced density matrix. One obtains [l6| 



~ L ■ 1 




— sin 




7T ' 





Si, 



(4) 



where si is a constant. As may be seen in Fig. [3l the 
numerically extracted central charge of the MC phase 
yields c = 1, as one would expect for a free boson, with 
coexisting gapped degrees of freedom. In addition, the 
AC+MC phase also has c = 1 . Note that it is not c = 2 
as would be the case for two independent Luttinger liq- 
uids. This reflects the coupled nature of the atomic and 
molecular condensates in the AC+MC phase, with ad- 
ditional gapped Ising degrees of freedom; the Feshbach 
term is relevant and drives the Z2 sector massive. Close 
to the MC to AC+MC transition, where the anticipated 
Ising gap closes, one expects the central charge to in- 
crease to c = 3/2, due to additional critical Ising degrees 
of freedom with c = 1/2. This is confirmed by DMRG 
in Fi g. [3] Further evidence is obtained from the differ- 
ence [lj AS(L) = S L (L/2) - S L/2 (L/4) = f ln(2) + . . . , 
as a function of e m ; see Fig. [3J For a given L this dis- 
plays a peak, whose location coincides with the MC to 
AC+MC transition obtained via the single-particle gap 
in Fig. [T] The evolution with increasing L is consistent 
with the passage towards c — 1 in the supcrfluid phases, 
and c = 3/2 in the vicinity of the transition. Applica- 
tion of this method to the MI to supcrfluid transitions 
Q yields c = 1 close to the MI boundary, suggesting XY 
behavior; see inset of Fig. |3l The absence of criticality 
within the MI phase is evidence against a super-Mott 
state and correlations decay exponentially [H, [l8| . 

Having provided evidence for a Z2 superfluid tran- 
sition, we now extract the Ising correlation length, £, 
and order parameter, ((f)), via finite size scaling of the 
atomic and molecular correlations. Due to the absence 
of particle conservation, and the presence of additional 
superfluid degrees of freedom, these cannot be readily 
obtained from the energy spectra alone. Ising scaling 
close to the transition implies that £ _1 ~ \M — M c \ and 
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FIG. 3. Top: Entanglement entropy Sl(1) obtained by 
DMRG in a periodic system with L — 64. We transit through 
the AC to AC+MC transition in Fig.[TJwith U = 0.7. The fits 
to Eq. Q yield c « 1 in the MC and AC+MC phases, and 
c f» 3/2 close to the transition. This reflects additional crit- 
ical Z,2 degrees of freedom. Due to the asymptotic nature of 
Eq. (T4J), high quality fits are obtained from the central region 
away from the boundaries. Bottom: Entanglement entropy 
difference AS(L) showing an Ising transition at e m ~ —3.8 
for U = 0.7, in agreement with Fig. [1] The solid lines are 
spline fits. Inset: A5* on passing through the AC+MC to MI 
transition at e m = 4, suggesting an XY transition with c = 1. 



((f)) ~ \Jvl — A^cl 1 / 8 , where M. is a mass scale parametriz- 
ing the departure from criticality. We identify the molec- 
ular density, A4 ~ p m , as the appropriate scaling vari- 
able. As shown in Fig. [4j the DMRG results arc in ex- 
cellent agreement with Ising critical exponents. This is 
non-trivial since the Ising degrees of freedom are non- 
local with respect to the atoms and molecules themselves. 
In summary, we have studied bosons interacting via Fes- 
hbach interactions in a ID lattice. We provide evidence 
for an Ising quantum phase transition between distinct 
supcrfluids. We extract both the Z2 order parameter, 
((f)), and the Ising correlation length, £. It would be in- 
teresting to see if this Z2 transition may be driven first 
order, and the effect of higher bands [19(- One may also 
consider Zn transitions involving N-particle pairing. 
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FIG. 4. DMRG with open boundaries. We use the parameters 
in Fig. [T] with U = 0.7 and transit from MC to AC+MC. 
Z2 disordered MC phase: (a) £ extracted from the atomic 
correlations with L = 256. £ _1 ~ \Pm— Pm\ v , where p m (pm ~ 
0.85) is the (critical) density of molecules and v = 1. Z2 
ordered AC+MC phase: (b) (</>) 2 , up to a constant prefactor, 
obtained by finite size scaling of the atomic correlations with 
up to L = 512. (c) Re-plotting yields /3 = 1/8. 
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